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The negatively charged nitrogen-vacancy (N-V~) color center in diamond is an important solid-state single 
photon source for applications to quantum communication and distributed quantum computation. Its fuU 
usefulness relies on sufficient radiative emission of the optical photons which requires realizable control to 
enhance emission into the zero-phonon line (ZPL) but untU now is still a challenge. Detailed understanding 
of the associated excitation process would be of essential importance for such objective. Here we report a 
theoretical work that probes the spin-conserving optical excitation of the N-V~ center. Using 
density-functional-theory (DFT) calculations, we find that the ZPL and the phonon-side band (PSB) depend 
sensitively on the axial strain of the system. Besides, we find a relatively small PSB appearing at about 
100 GPa in the emission spectrum at low temperatures, which provides a means to enhance the coherent 
emission of the N-V~ center in quantum optical networks. 

The coherent coupling between quantum objects and optical photons is of crucial importance in quantum 
information science' The associated quantum entangled states are widely used for applications to 
quantum information processings'^ ". Among the solid-state hosts which are preferred for these applica- 
tions because of scalability reasons, the N-V" center in diamond is a realizable solid-state quantum logic bit 
(qubit) as it offers long spin coherence time'"* '^, and the spin states can be optically initialized and read out""". 
Due to high photostability the N-V" center also presents a good candidate as a single-photon source used in 
quantum networks^". However, successful realization of long-distance spin-photon entanglement is still a chal- 
lenge. One of the difficulties arises form the ZPL emission of this defect including very little part of the total 
emission^'; while the frequency shifted PSB normally gives rise to deterioration of the spin-photon entanglement^. 

The N-V center in diamond consists of a N atom and a first-neighbor vacancy in the C site (shown in Fig. la). 
The electron paramagnetic resonance (EPR) experiments have revealed that the ground states and excited states 
involved in the optical excitations are both spin triplet (S = l)^^'^"*, implying an even number of electrons. 
Therefore, this defect is normally considered to be negatively charged and the extra electron is donated from 
the isolated N atom^*. For the negatively charged N-V" center with trigonal (Csv) symmetry, there are two fully 
symmetric fli states (labeled u and v) localized about the valance band maximum (VBM) and doubly degenerate e 
states (labeled fix and e^) deeply in the gap, with the occupation of six electrons^^: 
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Figure 1 | Structure of the N-V center in diamond, (a), Atomic structure 
of the N-V" center. The vacancy is indicated by a small yellow ball and the 
adjacent C and N by gray and blue balls, respectively (b), Electronic 
structure of the N-V" center: the calculated spin-polarized single- electron 
levels (with respect to the VBM), in ^Aj ground states and excited states. 

Where CTi_3 and CT4 are sp^ dangling bonds of the first-neighbor C and 
N atoms around the vacancy, respectively. ?^ is a parameter that 
determines the extent to which 04 is mixed in the v states, S is the 
overlap integral. 

In the ground states, the Ui states are fully occupied and the e states 
are occupied by two electrons of parallel spin, as illustrated in Fig. lb. 
Under the point group of C^y, the total wave function of the ground 
states has 'A2 symmetry. Upon optical excitation, as one electron in 
the spin-down v level is excited to the 'e^ or Cy level, the total spin of 
the system is conserved (shown in Fig. lb), and the resulting excited 
states is a spin triplet with the wavefunction of symmetry. 
Importantly, the change of the electronic configuration is followed 
by ionic relaxations, which affect strongly the optical features as the 
relaxation of atoms varies the transition energy, leading to large PSB 
both in absorption and emission. The large PSB in optical spectrum 
of this defect has been observed since few decades ago^', but the 
theoretical understanding of the associated mechanism remains a 
problem^'', as its interpretation is complicated by different results 
obtained for absorption and emission. 

It is possible for there to be strain in diamond, namely, the axial 
strain and non-axial strain. The latter lowers the symmetry of the N- 
V~ center, splits the excited states according to their orbital wave- 
functions and affects consequently the selection rules for the optical 
transitions'"'^'"^'-'. We don't discuss more here for the non-axial strain 
as it is not physically applicable to our objective. For the axial strain, it 
retains the N-V" center's Ct,v symmetry and is expected to shift the 
energy levels both in ground and excited states, varying the optical 
transition energies^". Furthermore, it is also possible to affect the 



oscillator strength of ZPL transitions, because of the strain changing 
the lattice parameter and lattice relaxation in the optical excitation. 
Until now, however, a complete picture for the dependence of the 
spin-conserving excitation on the axial strain of the lattice is still 
lacking. A detailed study to understand the excitation processes 
associated with the axial strain and lattice relaxation is therefore 
needed and is to be performed in this paper. 

Our theoretical work was carried out using DFT, which is a stand- 
ard method to investigate the properties of defect in solids. The 
geometry of the N-V center is modeled within a 216-atom cubic 
supercell of diamond. The axial strain of lattice was introduced by 
placing the N-V" center superceU in a hydrostatic pressure filed, 
while the strain was directly tuned through varying the pressure. 
For the excited states calculation, a preferred approach is to solve 
the Bethe-Salpeter Equation (BSE)"", together with employing a GW 
approximation^'. However, this approach is computationally prohib- 
itive for the calculations of geometric structures, and of electronic 
structures of large supercells. A realizable solution to this problem is 
the application of constrained DFT (CDFT)'"', which has been very 
successful in describing the excited states of solid-state system'^''* 
and was used here as the method to deal with the excited states of 
the N-V" center. Within CDFT, the calculated ZPL for the N-V" 
center in the absence of external strain is about 1.658 eV, which is 
14.7% smaller than the experimentally measured 1.945 eV^' but con- 
sistent with the previous CDFT-LDA'^' or CDFT-PBE results'^. This 
can be understood as a common self- interaction error of DFT"". 
Recently, Gali et al."" and Weber et al.""' have reported that using a 
screened Hartree-Fock hybrid density functional, more accurate ZPL 
for N- V~ center can be obtained. However, our goal in this paper was 
to give approximately the excitation energies to predict and under- 
stand the dependence of the optical excitation on the axial strain of 
lattice. We chose here the traditional density functional since it has 
been very successful in predicting the electronic'^ and vibration''' 
properties of N-V" center in diamond, and the pressure dependences 
for excited-states-related properties of solids, such as pressure coef- 
ficients of band gaps"* and pressure-induced metallization''^ etc. 

By minimizing the total energy of the N-V" system as a function of 
the ionic coordinates, we obtained the configuration coordinates {q^ 
for 'A2 ground states, then we placed the supercell in the stress field, 
and within the same way we obtained ground states of the 
N-V" center with the lattice strain under high pressures. Based on 
CDFT, we set the orbital occupancies for '£ states and fixed the 
electronic configuration while relaxing the ionic structure for each 
system obtained above to find the configuration coordinates for '£ 
excited states (q^), and then we calculated the potential energy curves 
(PEC) of the '£ states and 'A2 states along the ionic relaxation paths, 
as shown in Fig. 2a. We note that the PEC of the '£ states becomes 
constantly steeper as the external pressure increases, while the PEC of 
'A2 states flattens incipiently and then becomes steeper, revealing 
different changes of the potential energy surface (PES). As we wiU 
see, such difference gives rise to complicated changes of the excita- 
tion spectrum. 

Within the Franck-Condon approximation, the excitation spec- 
trum of the N-V" center can be described as follows: assuming that 
the transition of the electron is very fast compared with the relaxation 
of ions in the lattice, the transitions A ^ B and C ^ D (see Fig. 2a) 
are processes corresponding to the vertical absorption and emission 
where the lattice structures do not change. After the lattice relaxation 
taking place with the participation of phonons, the transitions 
between the energy minima of the 'A2 states and '£ states (A<-^C) 
can be realized, with a ZPL both in absorption and emission. The 
energy difference between the vertical absorption/emission and the 
ZPL corresponds to the relaxation energy, namely the Stokes-shift/ 
anti-Stokes shift (SS/ASS), which in experiment associates with the 
frequency shifted PSB of the highest intensity in absorption/emission 
at low temperatures. The calculated SS and ASS at zero pressure are 
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Figure 2 | Potential energy curve and excitation spectrum of tlie N-V~ 
center under liigli pressures, (a), Potential energy curves of ttie N-V" 
center under zero and iiigli pressures in the ground states and ^£ exited 
states along the ionic relaxation path in the states, and correspond 
to the energy minima of the ground and excited states, respectively. 

(b) , Calculated transition energies: vertical absorption and emission, 
and ZPL energies of the N-V" center under zero and high pressures. 

(c) , Calculated relaxation energies: Stokes Shift and Anti-Stokes Shift of the 
N-V" center under zero and high pressures. 

0.206 and 0.177 eV, respectively, comparable with the 0.235 and 
0.185 eV measured from experiment^'. The calculated relaxation 
energies in the optical excitation are considered to be more accurate 



than the transition energies, as they are obtained as the energy dif- 
ference between electronically similar conditions, and therefore the 
self- interaction errors are expected to be smaller. The obtained trans- 
ition energies, SS and ASS under external high pressures are shown in 
Fig. 2b and Fig. 2c, respectively. 

The results show very clearly that the vertical absorption energy 
rises significantly as the external pressure increases, while the vertical 
emission energy versus pressure curve rises incipiently and then 
becomes flattened and decays at ultra-high pressure (about 
300 GPa); meanwhile, the increasing of pressure also raises the 
energy of SS, but lowers the ASS in the range of 0-100 GPa and then 
raises it in higher range. The net result of these changes gives an 
increased ZPL. The obtained pressure coefficient for the shift of 
ZPL is about 5.75 meV/Gpa, in good agreement with the 5.5 meV/ 
Gpa measured from experimenf" (by Kobayashi et al.), in which the 
photoluminescence (PL) spectrum of the N- V" center was studied in 
a small pressure range (0-8.8 GPa). One difference between their 
experiment and our theoretical study is that they found pressure- 
independent PSBs, which significantly contradict our results. The 
reason for this is that the ZPL and PSB data summarized in figure 2 
of Ref 40 are not consistent with the PL spectrums given in the same 
work (in figure 1). As a result, they failed to recognize the change of 
the PSB due to pressure. On the contrary of their claim, we have 
checked the change of the PL spectrum in figure 1 of their work and 
then concluded that, e.g., the peak separation between the ZPL and 
the first phonon line decreases by 8% and 11% when the pressure 
goes from 0 to 5.1 and 8.8 Gpa, and this is in agreement with our 
calculated 9% and 15% decline of the ASS when going from 0 to 5 and 
10 GPa, respectively. Furthermore, the calculated ASS also agrees 
well with the decreased Huang-Phys factor obtained by Kobayashi 
et al. in the same work, revealing relatively suppressed PSB in this 
pressure range. 

To further understand the vertical absorption energies obtained, 
we calculated the spin-polarized single-electron levels of the N-V" 
defect in the states (see Fig. 3a). With the lattice strain under high 
pressure, the u level which appears in the valence band, together with 
the V level localized in the gap, shift down in energy with respect to 
the valence band maximum (VBM), lowering the spin-splitting ener- 
gies of the both. Simultaneously, the doubly degenerated e and e 
levels wUl shift up, resulting in larger energy range between the 
excitation-related levels. 

The defect levels depend not only upon the ionic environment but 
the electronic configuration as well. Therefore, either the change of 
the lattice arrangement or the varying of the orbital occupancies will 
give correspondingly orbital relaxations. This can be well verified by 
the calculated spin-polarized density of states (DOS, see Fig. 4) assoc- 
iated with the excitation: In the process of vertical absorption (A 
B), as one electron is excited from the v level to the e level, the spin up 
e states shift up in energy, reducing the spin-splitting energy with 
respect to the e states, and the spin up v states shift down to a lower 
energy [comparable to the vertical absorption, inverse orbital relaxa- 
tions of V and e states can be observed in the process of vertical 
emission (C — > D)]. Subsequently, the ionic relaxation (B C) will 
lower the e and e states, while raise the v and v states (we make no 
discussion for the u states as they are not sensitive to the ionic 
relaxation here). Because the v level is unoccupied, the e level is 
occupied and the e level is partially occupied in the ^£ states, such 
orbital relaxation consequently stabilizes the excited states of N- V" 
system; but the case is different for the deexcitation (D A) as the 
ionic relaxation takes place in the electronic configuration of ^A-i 
states, in which the orbital occupancies differ from the states 
as indicated in Fig. lb: the v level is occupied but the e level is 
unoccupied. Therefore, the energy difference associated with the 
orbital relaxation in the deexcitation would mostly cancel out, and 
the obtained ASS is supposed to originate mostly from the incidental 
deformation of the delocalized states created from the back bonds of 
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Figure 3 | Defect levels and atomic configuration under external high 
pressures, (a), Energy versus pressure curve for spin-polarized single- 
electron levels (with respect to the VBM) of the N-V" center in the 'A2 
states, (b), Calculated C-C bond length (Bc-c) in diamond; C-V (Dc-v) 
and N-V (Dn_v) distances of the N-V" center under high pressures in 'A2 
and states. 

the center, which are delocalized in valence and conduction bands of 
the host. This interprets the observed difference of the PEC for 
and ^£ states under high pressures. 

Once created in the bulk diamond, the N-V" center atoms will 
relax outward from the vacant site, the resulting extension strain 
retains the N-V" center's C^v symmetry and keeps the dangling 
bonds around the vacancy being highly localized. The calculated 
carbon-vacancy (C-V) and nitrogen -vacancy (N-V) distances are 
1.645 A and 1.698 A compared with the 1.550 A of the carbon- 
carbon (C-C) bond length in diamond (see Fig. 3b). Upon excita- 
tion, the C atoms further relax outward from the vacant site and 
the location of the N atom becomes less outward; the calculated 
displacement of C here is about 0.060 A, which is comparable with 
that of N (0.061 A). Because the defect center consists of three C 
but one N atom, the excitation will further lower the symmetry of 
the environmental tetrahedron (T^) crystal field and give even 
larger deformation of the delocalized states, and this way the 
ASS is created. We note that under high pressures, the displace- 
ments of the C atoms remain almost constant in the excitation, 
while the displacement of the N atom gets actually larger. Due to 
the lattice shrinkage under high pressure, in the reciprocal space, 
the defect levels will become more sensitive to the excitation. Fig. 4 
demonstrates that the orbital relaxations between the ground states 
and excited states (B C) become more evident under higher 
pressures, which is responsible for the steeper PEC of the states 
and the increasing SS obtained. 



As the pressure rises, the C-V and N-V distances decrease more 
rapidly than the C-C bond length, as shown in Fig. 3b. This will lower 
the local strain of the center. For instance, at 100 GPa, the calculated 
C-V distances are 1.464 A, which happens to be equal to the C-C 
bond length revealing almost no strain of the C, while the calculated 
N-V distance is 1.546 A revealing also a suppressed outward strain of 
N (0.082 A). This is because the interactions between the center's sp^ 
dangling bonds are not as strong as the bonding of the C atoms in 
diamond. Furthermore, the pressure will lower the relative strain 
between the ground and excited states if is not sufficiently high to 
induce an inward strain, lowering the ASS. This interprets the flat- 
tened PEC of the ^2 states and the decreasing ASS in the pressure 
range of 0-100 GPa. 

Note also that at 100 GPa, the calculated ASS converges very close 
to zero — — 1 3 me V. The obtained negative value is due to mosfiy the 
static J-T distortion in the excited states, as the geometric structure of 
the N-V~ system is fully relaxed without symmetry fixing. However, 
in most experimental cases (even at low temperatures), the breaking 
of the C3V symmetry due to J-T effect can not be observed"", as the 
potential well associated with the J-T distortion is of magnitude 
smaller than the zero-point energy^*, and the minima of the potential 
surface are normally occupied by the zero-point vibrations or higher 
vibrational states. Thus, from this point of view, the negative value 
obtained here also rises from the zero-point vibration (ZPV) error 
within the adiabatic approximation in DFT. To exclude the assoc- 
iated error, here we have tested the ASS in the absence of J-T distor- 
tion, the obtained value is 2 meV which again reveals a small ASS. 
The reason that the N-V" center has significantly small ASS at this 
pressure can be understood as follows: In the excited states, as the 
ionic relaxation takes place, the C-V and N-V distances shift to about 
1.525 and 1.461 A (see Fig. 3b), corresponding to outward strain of C 
(0.061 A) and very small inward strain of N (0.003 A), respectively. 
That is, after excitation, the strain of the N atom is correspondingly 
removed accompanied with the presence of smaller strain of the 
three C atoms, which tends to cancel out the energy difference 
between and in the 'A2 states and thus lowers the ASS. This 
can be understood also as a result of redistribution of the localized 
phonon states, which lowers the oscillator strength of the phonon 
emission and improves emission into the ZPL. Therefore, one would 
expect the N-V" center to display a relatively small emission PSB 
under about this pressure at low temperatures. Although can not be 
avoided here, due to the flattened PEC of the 'A2 states, the intrinsic 
ZPV error is expected to be even smaller compared with the one at 
zero pressure, and this also holds for a small PSB of the emission 
spectrum. 

The inward strain of the center wiU take place and become evident 
as the pressure goes up further (see Fig. 3b) and consequently 
increases the relative strain between the ground and excited states. 
The DOS obtained here reveal different deformations of the deloca- 
lized states, which are noted to be relatively significant at 500 GPa 
and evidence the steep PEC of the M2 states and the large ASS 
obtained here. Note that the v levels become somewhat delocalized 
along with the inward strain. The calculated DOS show clear broad- 
ening of the V states at 300 and 500 GPa. This results from close 
interactions between the center's atoms as the distances between 
the neighboring C become relatively shorter, and the overlap 
between the sp'' dangling bonds can be no more negligible. The v 
levels shift down in energy, couple strongly to the VBM states and 
become somewhat delocalized, lowering the energy. Indeed, we note 
that the delocalization of the v levels is not as significant in the excited 
states due to the ionic relaxation, which is also responsible for the 
large ASS and the inversely changed vertical emission energy starting 
at about 300 GPa. 

In fact, the relaxation energy in the optical absorption/emission of 
the N-V" center arises from the orbital relaxation of the defect 
together with the deformation of the delocalized states, which are 
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Figure 4 | The density of states (DOS) of the N-V~ center, (a), Calculated spin-polarized DOS of the N-V" center at 0 and 100 GPa. 
(b), Spin-polarized DOS at 300 and 500 GPa. A, B, C, and D correspond to the ionic and electronic configurations indicated in Fig. 2a. The black 
and blue lines represent the spin-up and -down states; the vertical gray lines show the corresponding defect levels, where the solid, dot and dash lines 
indicate occupied, unoccupied and half occupied orbitals, respectively; the horizontal arrows represent the associated orbital shifts in the excitation/ 
deexcitation processes. The zero of energy is set at the Fermi level of the ground states. 



both correlated strongly with the local ionic arrangement. In the 
optical excitation, the relaxation energy associated with the orbital 
relaxation is relatively larger than that from the deformation of delo- 
calized states and thus dominates the SS; whereas in the deexcitation, 
due to different electronic occupancies the ASS is derived mostly 
from the latter. This provides an understanding of how ionic relaxa- 
tion affects the excitation/deexcitation of the N- V" system and inter- 
prets the experimental observation^'. The obtained results also 
demonstrate a dependence of the excitation spectrum of the N-V" 
center on the axial strain of the lattice, which is directly tunable via 
the external pressure. We expect based on these results that pressure- 
control for the optical transitions of the N-V" center will be possible. 
Furthermore, the relatively small ASS obtained at 100 GPa demon- 
strates a suppressed PSB of the emission spectrum at low tempera- 
tures. This might prove useful as a means to improve the ZPL 
emission of the N-V" center for applications to quantum cryp- 
tography and quantum communication. Beyond these specific appli- 
cations, the realizable control of the PSB emission also paves a 
pathway towards improving the efficiency of photon-spin coherent 
coupling, which is currently the key limitation in realization of long- 
distance photon-spin entanglement. 

Methods 

We optimized the geometry and calculated the charge density for the supercell 
structure of N-V" center employing spin polarized generalized gradient approxi- 
mation (SGGA) with Perdrew-Burke-Enzerhof (PBE)''^ exchange-correlation 



functional. Projected augmented wave (PAW) pseudopotentials*"* for the C and N 
were used with a planewave basis-set cutoff energy of 520 eV, which is adequate for 
optimization of the geometry under high pressure. We used a 3 X 3 X 3 Monkhorst- 
Pack mesh for the k point samplings in the Brillouin zone. This provides a very well 
converged charge density. The 216-atom supercell of diamond was constructed by 
extending the conventional cubic cell three times along the lattice vector a, b and c, 
respectively; we then placed the N-V defect in the center of the supercell to model the 
defective system. The optimized lattice constant of the supercell is about 10.73 A at 
0 Pa and 9.11 A at 500 GPa, which is appropriate for simulation of the isolated N-V 
defect and is also large enough to obtain accurately the associated spin density. Within 
CDFT, the electronic occupancies for the excited states were obtained by setting 
zero occupation for the spin-down v level and one occupation for spin-down or Cy 
level. In the geometric optimization, no symmetry restriction was applied, with all the 
atoms allowed to relax untH the calculated Hellmann-Feynman forces became smaller 
than 1 meV/A. 
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